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Abstract. A fully relativistic modified gravitational theory including a fifth force 
skew symmetric field is fitted to the Pioneer lO/II anomalous acceleration. The theory 
allows for a variation with distance scales of the gravitational constant G, the fifth 
force skew symmetric field coupling strength u; and the mass of the skew symmetric 
field /i = 1/A. A fit to the available anomalous acceleration data for the Pioneer 10/11 
spacecraft is obtained for a phenomenological representation of the "running" constants 
and values of the associated parameters are shown to exist that are consistent with 
fifth force experimental bounds. The fit to the acceleration data is consistent with all 
current satellite, laser ranging and observations for the inner planets. 
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1. Introduction 

The radio tracking data from tfie Pioneer 10/11 spacecraft during tlieir travel to tlie 
outer parts of the solar system have revealed an anomalous acceleration. The Doppler 
data obtained at distances r from the Sun between 20 and 70 astronomical units 
(AU) showed the anomaly as a deviation from Newton's and Einstein's gravitational 
theories. The anomaly is observed in the Doppler residuals data, as the differences of the 
observed Doppler velocity from the modelled Doppler velocity, and can be represented 
as an anomalous acceleration directed towards the Sun, with an approximately constant 
amplitude over the range of distance, 20 AU < r < 70 AU [1, 2, 3]: 

ap = (8.74 ± 1.33) x 10"^ cm s'^. (1) 

After a determined attempt to account for all known sources of systematic errors, the 
conclusion has been reached that the anomalous acceleration towards the Sun could be 
a real physical effect that requires a physical explanation [1, 2, 3]. 

Two theories of gravity called the metric-skew-tensor gravity (MSTG) theory [4] 
and the scalar-tensor vector gravity (STVG) theory [5] have been proposed to explain 
the rotational velocity curves of galaxies, clusters of galaxies and cosmology. A fitting 
routine for galaxy rotation curves has been used to fit a large number of galaxy rotational 
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velocity curve data, including low surface brightness (LSB), high surface brightness 
(HSB) and dwarf galaxies with both photometric data and a two-parameter core model 
without non-baryonic dark matter [5, 6]. The fits to the data are remarkably good and 
for the photometric data only the one parameter, the stellar mass-to-light ratio, (M/L), 
is used for the fitting, once two parameters Mq and Vq are universally fixed for galaxies 
and dwarf galaxies. A large sample of X-ray mass profile cluster data has also been 
fitted [7]. 

The gravity theories require that Newton's constant G, the coupling constant u that 
measures the strength of the coupling of the skew field to matter and the mass fi of the 
skew field, vary with distance and time, so that agreement with the solar system and the 
binary pulsar PSR 1913-1-16 data can be achieved, as well as fits to galaxy rotation curve 
data and galaxy cluster data. In Ref. [5] and Ref. [8], the variation of these constants 
were based on a renormalization group (RG) flow description of quantum gravity theory 
formulated in terms of an effective classical action. Large infrared renormalization effects 
can cause the effective G, u, ^ and the cosmological constant A to run with momentum 
k and a cutoff procedure leads to a space and time varying G, uj and yU, where /i = l/rg 
and ro is the effective range of the skew symmetric field. In the STVG theory [5], 
the action contains contributions that lead to effective field equations that describe 
the variations of G, /i and the coupling constant uj that measures the strength of the 
coupling of the skew symmetric field B^^, with matter. In principle, we can solve for the 
complete set of field equations and determine the dynamical behavior of all the fields. 
However, in practice we make approximations allowing us to obtain partial solutions to 
the equations, yielding predictions for the various physical systems considered. 

Both the MSTG and STVG theories lead to the same modified acceleration law 
for weak gravitational fields and the same fits to galaxy rotation curve and galaxy 
cluster data, as well as to agreement with the solar system and pulsar PSR 1913+16 
observations. However, the STVG theory is simpler in its structure, so in the following, 
we shall restrict our attention to this theory. 

An important constraint on gravity theories is the bounds obtained from weak 
equivalence principle tests and the existence of a "fifth" force, due to the exchange of 
a massive vector boson [9, 10, 11, 12, 13, 14]. These bounds cannot rule out gravity 
theories that violate the weak equivalence principle or contain a fifth force at galactic 
and cosmological distance scales. However, we shall in the following attempt to explain 
the anomalous acceleration observed in the Pioneer 10/11 spacecraft data [1, 2, 3] and 
simultaneously fit the solar planetary data. For our study of the solar system, we must 
account in our modified gravity theory not only for the variation of G with respect to 
the radial distance r from the center of the Sun, but also for the variations of /i = l/rg 
and UJ with respect to r. Since we do not possess rigorous solutions for the variations 
of ro and uj, we use a phenomenological parameterization of the varying parameters 
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tq and 00 to obtain fits to the anomalous Pioneer acceleration data that are consistent 
with the solar system and fifth force experimental bounds. As has been demonstrated 
by Ref. [15], the fifth force experimental bounds rule out a gravitational and fifth force 
explanation of the Pioneer anomaly for fixed universal values of the coupling strength 
a and the range A = tq. 

In an asymptotically free RG fiow description of quantum gravity, a distance 
(momentum) scaling law for the gravitational coupling constant G, the coupling constant 
u and the range A leads to strong infrared renormalization effects for increasing distance 
scale [8, 4]. We shall find in the following that an asymptotically free-type of gravity 
and fifth force distance scaling behavior can possibly explain the anomalous Pioneer 
acceleration and be in agreement with local Eotvos experiments and planetary deep 
space probes. 

An important feature of the MSTG and STVG theories is that the modified 
acceleration law for weak gravitational fields has a repulsive Yukawa force added to 
the Newtonian acceleration law for equal "charges" . This corresponds to the exchange 
of a massive spin-1 boson, whose effective mass and coupling to matter can vary with 
distance scale. A scalar component added to the Newtonian force law for equal charges 
would correspond to an attractive Yukawa force and the exchange of a spin-0 particle. 
The latter acceleration law cannot lead to a satisfactory fit to galaxy rotation curves 
and galaxy cluster data [6, 7]. 

All the current applications of the two gravity theories that can be directly 
confronted with experiment are based on weak gravitational fields. To distinguish the 
MSTG and STVG theories, it will be necessary to obtain experimental data for strong 
gravitational fields, for example, black holes as well as cosmology. 

2. Weak Fields and the Modified Gravitational Acceleration 

The equation of motion of a test particle, about a mass M, obtained for weak 
gravitational fields and for a static spherically symmetric solution of the massive skew 
field Bfj^i, in the STVG theory takes the form [5]: 

+ + (2) 

where K and fi are positive constants and Jn is the Newtonian angular momentum. We 
have not assumed in Equation (2) any particular composition dependent model for K 
for the additional Yukawa force term. 

We observe that the additional Yukawa force term in Equation (2) is repulsive for 
like charges in accordance with the exchange of a spin-1 ~ massive boson. We shall treat 
a = a{r) and fi = 1/A(r) as running parameters which are effective functions of r and 
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obtain for the radial acceleration on a test particle 

.M^-%^ + i^(r)!2^<^(l + 3^). (3) 

Here, we are treating the "constants" K and A as running with distance scale, in accord 
with our effective RG flow interpretation of G(r), K{r) and A(r) [4, 8]. On the other 
hand, in the STVG theory [5] the varying constants are treated as scalar fields and are 
solved in terms of effective field equations obtained from an action. In Equation (3), 
Goo denotes the renormalized value of the gravitational constant: 



Goo = Gq{i + ac 



where Go = 6.674 x 10 ^g ^cm^s ^ is Newton's "bare" gravitational constant. Let 
K{r) = GoMa{r). 
By using Equation (4), we can rewrite the acceleration in the form 
, , GoM 



1 + a(r) [l - exp(-r/A(r)) f 1 + 

I V A(r)/ 



The acceleration law (6) can be written 
G(r)M 



a{r) 



where 



G(r) = Go [l + air) (l - exp(-r/A(r)) (l + ^) ) 
is an effective expression for the variation of G with respect to r. 



(4) 
(5) 
(6) 

(7) 
(8) 



3. Pioneer Anomalous Acceleration 



We postulate a gravitational solution that the Pioneer 10/11 anomaly is caused by 
the difference between the running G(r) of Equation (8) and the bare value. Go = 
6.674 X lO^^g^^cm^s"^. So the Pioneer anomalous acceleration directed towards the 
center of the Sun is given by 



ap 



6G{r)MQ 



where 



^^(r) = Goa(r) 1 - exp(-r/A(r)) ( 1 + 



A(r) 



(9) 



(10) 



We propose the following parametric representations of the "running" of a(r) and A(r): 
«(r) = aoo(l-exp(-r/f))^/^ (11) 

Ano 



A(r) 



(1 — exp(— r/r))' 



(12) 
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Figure 1. Best fit to tlic Pioneer 10/11 anomalous acceleration data extracted from 
Figure 4 of Ref. [16] plotted against the position, r in AU, on a logarithmic scale out 
to r = 10,000AU. Pioneer 10 data is shown with open cyan circles and Pioneer 11 
data is shown with closed blue circles. The Pioneer 10/11 anomalous acceleration, ap, 
is in units of 10^^ cm/s and the constant Pioneer anomaly result of Equation (1) is 
shown with dotted black lines. 



Here, f is a non- running distance scale parameter and 6 is a constant. 

In Figures 1 and 2, we display a best fit to the acceleration data extracted from 
Figure 4 of Ref. [16] obtained using a nonlinear least-squares fitting routine including 
estimated errors from the Doppler shift observations [2]. The best fit parameters are: 

Ooo = (1.00 ±0.02) X 10"^ 
Aoo = 47 ± 1 AU, 
f =4.6±0.2AU, 
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b =4.0. (13) 

The small uncertainties in the best fit parameters are due to the remarkably low variance 
of residuals corresponding to a reduced degree of freedom of 0.42 signalling a 

good fit. 

11 I I I I I 




r [AU] 



Figure 2. Best fit to the Pioneer 10/11 anomalous acceleration data extracted from 
Figure 4 of Ref. [16] plotted against the position, r in AU, on a linear scale out to 
r = 50 AU. Pioneer 10 data is shown with open cyan circles and Pioneer 11 data is 
shown with closed blue circles. The Pioneer 10/11 anomalous acceleration, ap, is in 
units of 10~^ cm/s and the constant Pioneer anomaly result of Equation (1) is shown 
with dotted black lines. 



Fifth force experimental bounds plotted for logj^o versus log^Q A are shown in 
Fig. 1 of Ref. [15] for fixed values of a and A. The updated 2003 observational 
data for the bounds obtained from the planetary ephemerides is extrapolated to 
r = 10^^ m = 6,685AU [17]. However, this extrapolation is based on using fixed 
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universal values for the parameters a and A. Since known reliable data from the 
ephemerides of the outer planets ends with the data for Pluto at a distance from 
the Sun, r = 39.52 AU = 5.91 x 10^^ m, we could claim that for our range of values 
47 AU < A(r) < oo, we predict air) and A(r) values consistent with the un- extrapolated 
fifth force bounds. 

A consequence of a variation of G and GMq for the solar system is a modification 
of Kepler's third law: 

alL = GiapL)Mj'^)\ (14) 



271 

where Tpi is the planetary sidereal orbital period and apL is the physically measured 
semi-major axis of the planetary orbit. For given values of apL and Tpi, Equation (14) 
can be used to determine G{r)MQ. The standard method is to use astrometric data to 
define GMq for a constant value, 

G(r)MQ = GMMq = k\ (15) 

where a® is the semi-major axis for Earth's orbit about the Sun, and n is the Gaussian 
gravitational constant given byf 

K = 0.01720209895 AU^^Vday. (16) 

We obtain the standard semi-major axis value at 1 AU: 

4^ = G(ae)MQ(^)'. (17) 

For several planets such as Mercury, Venus, Mars and Jupiter there are planetary ranging 
data, spacecraft tracking data and radiotechnical fiyby observations available, and it is 
possible to measure api directly. For a distance varying GMq we derive [9, 10]: 

-G(api)Mo] 1/3 



rapL\ 
Kapi) 



1 + r]PL 



1^ 



1. (19) 



Here, it is assumed that GMq varies with distance such that rjpi can be treated as 
a constant for the orbit of a planet. We may substitute the Gaussian gravitational 
constant of Equation (15) into Equation (18) and obtain 

\G{apL)y" 
Vpl - -jT, — T 

In Figure 3, we display the variation of 5G/Gq arising from Equation (10) versus 
r for the parametric values of a{r) and A(r) of Equation (11) and Equation (12), 
respectively, using the best fit values for the parameters given in Equation (13). The 
behaviour of G{r)/Go is closely constrained to unity over the inner planets until beyond 
the orbit of Saturn (r > 10 AU) where the deviation in Newton's constant increases to an 
asymptotic value of Goo/ Go 1.001 over a distance of hundreds of AU. The transition 

t http:/ /ssd.jpl. nasa.gov/astro_constants. html 
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Figure 3. Plot of 5G{r)/GQ arising from Equation (10) versus r in AU predicted by 
the best fit parameters of Equation (13) to the Pioneer 10/11 anomalous acceleration 
data. 

to the renormalized Newtonian constant is somewhat sudden as though the underlying 
physics was mimicking some source of dark matter. However, the fit is dictated by 
observations within the solar system and it is the nature of the renormalization group 
flow to predict a continuous phase transition in the behaviour of the coupling constant 
in the vicinity of the infrared fixed point. It is this underlying physics which we have 
applied phenomenologically to solve the problem of galaxy rotation curves without dark 
matter [6] and galaxy cluster masses without dark matter [7]. 

For the nine planets, we obtain the values of rjpL shown in Table 1. We see that 
we are able to obtain agreement well within the bounds of possible variation of GMq 
consistent with the data [9, 10] for Mercury, Venus, Mars and Jupiter. No observational 
limit on rjpL for Saturn or the outer planets has yet been established; but this is 
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Table 1. Theoretical predictions of the values of i]pl of Equation (19) compared to 
the observational limits taken from Ref. [10] and the best-fit theoretical predictions 
for the Pioneer Anomaly, ap, for the nine planets at a distance r from the Sun. No 
observational limits were computed beyond Saturn in Ref. [10] due to uncertainty in 





optical data. Beyond the outer planets, the theoretical predictions for T]{r) a 




the asymptotic value 7700 = 


3.34 X 10^**. 




Planet 
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iheoretical rredict 


ion (Jbservationai Limit 


ap 




(AU) 


VPL (10-1°) 


VPL (10 


(10 8 cm/s ) 


Mercury 


0.38 


-6.55 X 10-5 


+40 ± 50 


1.41 X 10-1° 


Venus 


0.72 


-6.44 X 10-5 


-55 ± 35 


5.82 X 10-8 


Earth 


1.00 


0.00 X 10° 





1.16 X 10-6 


Mars 


1.52 


4.93 X 10-3 


-0.9 ±2.1 


4.42 X 10-5 


Jupiter 


5.20 


4.19 X 10^ 


+200 ± 400 


2.76 X 10-1 


Saturn 


9.54 


1.67 X 10^ 




3.27 X 10° 


Uranus 


19.22 


1.84 X 10^ 




8.86 X 10° 


Neptune 


30.06 


4.39 X 10^ 




8.65 X 10° 


Pluto 


39.52 


6.77 X 10^ 




7.72 X 10° 



precisely where the deviation 5G{r)/GQ leads to a sizable contribution in the theoretical 
prediction for Vpl- The reason for the uncertainty beyond the orbit of Saturn and the 
lack of observational limits on rjpi is that the ephemerides for the outer planets is based 
on optical measurements. Even in the context of Newton's theory, the extrapolation 
of Kepler's third law of Equation (14) using the Gaussian gravitational constant of 
Equation (16) which fits the inner planets misestimates the semi-major axis, apL, or 
the orbital period, Tp^, of the outer planets resulting from Newtonian perturbations 
due to Jupiter and the gas giants and their satellites, the Kuiper belt and hundreds of 
asteroids. The latest version of the planetary part of the numerical ephemerides is a 
numerical integration of the post-Newtonian metric. It attempts to account for these 
perturbations from Kepler's law beyond Saturn by a least squares adjustment to all 
the available observations including the CCD optical astrometric observations of the 
outer planets. These values (without uncertainty) are available from the Solar System 
Dynamics Group (SSD) of the Jet Propulsion Laboratory (JPL) through the Horizon's 
ephermeris DE410 online|. The Russian Academy of Sciences has also placed their 
latest values known as EPM2004 online§. Because the perturbations change daily due 
to the motion within the solar system, the planetary ephemerides quoted values for api 
and TpL change daily. In order to compute deviations from Kepler's third law for the 

I http:/ /ssd.jpl. nasa.gov/horizons. html 

§ ftp://quasar.ipa.nw.ru/incoming/EPM2004 
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outer planets, we have listed today's best known values in Table 2. The uncertainty 
in the EPM2004 deduced values for the semi-major axes of the planets, Aap^, have 
been studied in Ref. [18] and the quoted values are listed in Table 2. Ref. [18] warns 
that the real errors may be larger by an order of magnitude. The uncertainty in the 
periods for the outer planets are not quoted in either EPM2004 or DE410, and so we 
have assumed small uncertainties based on the precision provided by the JPL Horizon's 
online ephemeris. We may calculate the uncertainty, Atjpl, by propagating the errors 
AapL and ATp^ according to Equations (17) and (18), neglecting any uncertainty in 
the Gaussian graviational constant of Equation (16): 

A,P. ^ J(^)\a^^^)\ (20, 

V V O-PL J \^apL TpL J 

Although according to Table 1 we are consistent with the observational limits 
of rjpL for theinner planets to Jupiter, the computation of Ref. [10] attempted to set 
model-independent constraints on the possible modifications of Newtonian gravity. The 
procedure was to run the planetary ephemerides numerical integration with the addition 
of TjpL as free parameters. Because there was one additional parameter for each planet, 
they were only able to find observational limits for the inner planets including Jupiter. 
In order to compute the observational limit for rjpL for the outer planets, it would be 
necessary to compute the planetary ephemerides using the modified acceleration law of 
Equations (7) and (8). Although this is beyond the scope of the current investigation, 
we may approximate here the observational limit of rjpi for the outer planets as the 



Table 2. Mean ephemerides orbital parameters for the semi-major axes (J2000), 
orbital eccentricities and the sidereal orbital periods (JPL Horizon's online ephemeris). 
The errors in the semi-major axes are deduced from Table 4 of Ref. [18] with 1 AU = 
(149597870696.0 ± 0.1) m. The computation of A?/pl is based on the propagation of 
the uncertainties according to Equation (20). 



Planet 


apL 


AapL 






epL 


TpL 


ATpL 




ArjpL 






(AU) 


(AU) 








(days) 


(days) 




(10-1° 


) 


Mercury 


0.38709893 


7.02 X 


10' 


-13 


0.206 


87.968435 


5.0 X 10" 


7 


3.79 X 


10^ 


Venus 


0.72333199 


2.20 X 


10- 


-12 


0.007 


224.695434 


5.0 X 10" 


7 


1.48 X 


10^ 


Earth 


1.00000011 


9.76 X 


10- 


-13 


0.017 


365.256363051 


5.0 X 10- 


10 


1.33 X 


10-2 


Mars 


1.52366231 


4.39 X 


10- 


-12 


0.093 


686.980 


5.0 X 10- 


4 


4.85 X 


10^ 


Jupiter 


5.20336301 


4.27 X 


10- 


-9 


0.048 


4330.595 


5.0 X 10" 


4 


7.68 X 


102 


Saturn 


9.53707032 


2.82 X 


10- 


-8 


0.056 


10746.94 


5.0 X 10" 


2 


3.10 X 


10^ 


Uranus 


19.19126393 


2.57 X 


10- 


-7 


0.047 


30685.4 


1.0 X 10° 




2.17 X 


10^ 


Neptune 


30.06896348 


3.20 X 


10- 


-6 


0.009 


60189. 


5.0 X 10° 




5.54 X 


10^ 


Pluto 


39.48168677 


2.32 X 


10- 


-5 


0.250 


90465. 


1.0 X 10^ 




7.37 X 


10^ 
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uncertainty A?7pi from Equation (20), for the perturbations of Figure 3, 6G{r)/GQ, are 
small compared to the Newtonian perturbations acting on the outer planets. The results 
for Ar]PL due to the uncertainty in the planetary ephemerides are presented in Table 2 
for the nine planets and exceed the predictions, rjPL, of Table 1. 

The relativistic equation of motion for a test particle in our gravitational theory 
may be solved perturbatively in a weak field approximation for the anomalous perihelion 
advance of a planetary orbit: 



where we have assumed as with Kepler's third law that GMq and a vary with distance 
such that they can be treated as constants for the orbit of a planet, where we have made 
use of the approximation G{r) ^ Gq [5], which is the case from the fit to the Pioneer 
10/11 anomalous acceleration data. We may rewrite Equation (21) as the perihelion 
advance in arcseconds per century: 

^PL= ^ rj. = -2 71 2 \rr , 22 

2'kTpl c^apiil - ef,jTpL 

where Tp^ is the planetary orbital period, and ep^, is the planetary orbital eccentricity. 
We may separate Equation (22) into the usual Einstein anomalous perihelion advance, 
and a prediction of the correction to the anomalous perihelion advance: 

ujpL = uJo + ^i, (23) 



Table 3. The values of the running parameters. a{r) of Equation (11) and A(r) of 
Equation (12) and the deviation in the dimensionless gravitational constant, 5G{r)/Go 
of Equation (8), calculated for each planet. Included on the right of the table is the 
theoretical (Einstein) perihelion advance of Equation (24), and the predicted retrograde 
of Equation (25) for the planets, and the limits set by the ephemeris [19]. 



Planet 


apL 




XpL (AU) 


SGpl/ Gq 


Einstein 


uj ("/century) 
Retrograde Ephemeris 


Mercury 


6.51 X 


10- 


-6 


1.11 X 10^ 


3.98 X 10^ 


19 


42.99 


-2.80 X 


10^ 


-4 


-0.0336 ±0.0050 


Venus 


2.12 X 


10- 


~5 


1.05 X 10^ 


5.04 X 10" 


16 


8.63 


-1.83 X 


10" 


-4 




Earth 


3.82 X 


10" 


-5 


3.23 X 10^ 


1.84 X 10" 


14 


3.84 


-1.47 X 


10- 


-4 


-0.0002 ± 0.0004 


Mars 


7.95 X 


10' 


-5 


7.44 X 10^ 


1.67 X 10^ 


12 


1.35 


-1.07 X 


10- 


-4 


0.0001 ±0.0005 


Jupiter 


4.59 X 


10' 


--4 


2.23 X 10^ 


1.23 X 10' 


-7 


0.0624 


-2.86 X 


10" 


-5 


0.0062 ±0.036 


Saturn 


7.64 X 


10" 


-4 


8.05 X 10^ 


4.96 X 10" 


-6 


0.0137 


-1.05 X 


10" 


-5 


-0.92 ±2.9 


Uranus 


9.69 X 


10- 


-4 


5.00 X 10^ 


5.55 X 10' 


-5 


0.00239 


-2.31 X 


10" 


-6 


0.57 ± 13.0 


Neptune 


9.97 X 


10" 


-4 


4.73 X 10^ 


1.34 X 10' 


-4 


0.00078 


-7.73 X 


10" 


-7 




Pluto 


1.00 X 


10' 


~3 


4.70 X 10^ 


2.05 X 10' 


-4 


0.00042 


-4.18 X 


10" 


-7 
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(24) 



ui = -apLUjQ 



(25) 



is the predicted retrograde (note the minus-sign in Equation (25)) to the anomalous 
perihelion advance of Equation (24). The measured perihelion precession is best known 
for the inner planets (for Mercury the precession obtained from ranging data is known 
to 0.5% [20]). For each of the planets in the solar system, we find that apL << 1, so 
that our fit to the Pioneer anomalous acceleration is in agreement with the relativistic 
precession data. The results for the Einstein perihelion advance, and our predicted 
retrograde for each planet, and the observational limits set by the recent ultra-high 
precision ephemeris are listed in Table 3. 

The validity of the bounds on a possible fifth force obtained from the ephemerides of 
the outer planets Uranus, Neptune and Pluto are critical in the exclusion of a parameter 
space for our fits to the Pioneer anomaly acceleration. Beyond the outer planets, the 
theoretical prediction for ri{r) approaches an asymptotic value: 



We see that the variations ("running") of a{r) and A(r) with distance play an important 
role in interpreting the data for the fifth force bounds. This is in contrast to the standard 
non-modified Yukawa correction to the Newtonian force law with fixed universal values 
of a and A and for the range of values < A < oo, for which the equivalence principle and 
lunar laser ranging and radar ranging data to planetary probes exclude the possibility 
of a gravitational and fifth force explanation for the Pioneer anomaly. 

4. Conclusions 

A modified gravity theory based an a. D = A pseudo-Riemannian metric, a spin-1 vector 
field and a corresponding skew symmetric field B^i, and dynamical scalar fields G, 
uj and yU, yields a static spherically symmetric gravitational field with an additional 
Yukawa potential and with effective varying coupling strength a{r) and distance range 
A(r). This modified acceleration law leads to remarkably good fits to a large number of 
galaxies [6] and galaxy clusters [7] without non-baryonic dark matter. The previously 
published gravitational theories MSTG [4] and STVG [5] yielded the same modified weak 
gravitational field acceleration law and, therefore, the same successful fits to galaxy and 
cluster data. The MSTG and STVG gravity theories can both be identified generically 
as metric-skew-tensor gravity theories, for they both describe gravity as a metric theory 
with an additional degree of freedom associated with a skew field coupling to matter. 



lim r]{r) = 3.34 x 10 



-4 



(26) 
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For MSTG, this is a third-rank skew field -F)jjyA, while for STVG the skew field is a 
second-rank tensor B^^. 

An action Ss for the scalar fields G{x), uj{x) and = 1/A(x) and the field 

equations resulting from a variation of the action, 6Ss = 0, can be incorporated into 
the MSTG and STVG theories. The dynamical solutions for the scalar fields give an 
effective description of the running of the constants in an RG fiow quantum gravity 
scenario [8, 4], in which strong infrared renormalization effects and increasing large 
scale spatial values of G and u lead, together with the modified acceleration law, to 
a satisfactory description of galaxy rotation curves and cluster dynamics without non- 
baryonic dark matter, as well as a solution to the Pioneer 10/11 anomalous acceleration 
data. 

We have demonstrated that the STVG theory can explain the Pioneer anomalous 
acceleration data and still be consistent with the accurate equivalence principle, lunar 
laser ranging and satellite data for the inner solar system as well as the outer solar 
system planets including Pluto at a distance of r = 39.52 AU = 5.91 x 10^^ meters. The 
ephemerides for the outer planets are not as well know as the ones for the inner planets 
due to their large distances from the Sun. The orbital data for Pluto only correspond 
to the planet having gone round 1/3 of its orbit. It is important that the distance range 
parameter lies in the region 47 AU < A(r) < oo for the best fit to the Pioneer acceleration 
data, for the range in the modified Yukawa correction to Newtonian gravity lies in a 
distance range beyond Pluto. Further investigation of fifth force bounds obtained by 
an analysis of the planetary data for the outer planets, based on the modified gravity 
theory is required. We are predicting that measurements of a fifth force in the solar 
system will become measurable at distances r > 10 AU from the Sun where as shown in 
Figure 3, 6G{r)/Go (and tipl) become potentially measureable. 

Perhaps, a future deep space probe can produce data that can check the predictions 
obtained for the Pioneer anomaly from our modified gravity theory. Or perhaps utilizing 
minor planets may clarify whether the Pioneer anomaly is caused by the gravitational 
field in the outer solar system [21]. An analysis of anomalous acceleration data obtained 
from earlier Doppler shift data retrieval will clarify in better detail the apparent onset of 
the anomalous acceleration beyond the position of Saturn's orbit. The planned LATOR 
mission will achieve significantly improved measurements of possible new degrees of 
freedom such as dynamical scalar and vector fields in the solar system, which could test 
the proposed MSTG and STVG theories [22]. 
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